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ABSTRACT 

We report the discovery of WASP-4b, a large transiting gas-giant planet with an orbital period 
of 1.34 days. This is the first planet to be discovered by the SuperWASP-South observatory 
and CORALIE collaboration and the first planet orbiting a star brighter than 16 th magnitude 
to be discovered in the Southern hemisphere. A simultaneous fit to high-quality lightcurves 
and precision radial- velocity measurements leads to a planetary mass of 1.22^0 oi Mj up an d a 
planetary radius of 1.42±°° 7 R Jup . The host star is USNO-B1.0 0479-0948995, a G7V star of 
visual magnitude 12.5. As a result of the short orbital period, the predicted surface temperature 
of the planet is 1761 K, making it an ideal candidate for detections of the secondary eclipse at 
infrared wavelengths. 



Subject headings: planetary systems: individual: WASP-4b 
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1. Introduction 

Transiting extrasolar planets offer our best 
opportunity for measuring fundamental param- 
eters of extrasolar planets, such as radius and 
mass, allowing us to test our models for plane- 
tary formation and evolution. To-date the only 
transiting planets known in the Southern sky are 
those disco vered by surveys targeting the galactic 
plane (e.g. 



Udalski et al 
20081) 



arg 
1 



2002; Pont et all 12007 



Weldrake et al.1 120081 ) and find planets around 
stars with typical visual magnitudes of 16-17. 

The SuperWASP-South (WASP-S) instrument 
looks for transiting planets around stars of visual 
magnitude 8-13, and will eventually cover the 
entire southern sky, excluding the Galactic plane. 
We report here the discovery of WASP-4b, a tran- 
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siting hot- Jupiter orbiting a star of magnitude 
12.5. 

This is the first discovery by WASP-S, made in 
collaboration with radial-velocity measurements 
from the CORALIE spectrograph on the 1.2-m 
Euler telescope. The discovery marks the begin- 
ning of a campaign to discover brighter transiting 
extrasolar planets in the Southern hemisphere, 
to complement those di scovered in the No rth by 
projects including HAT dBakos et alj|2002l), TrES 
(lO'Donovan et al.l | 2006h. XO dMcCullough et all 



20051 ) and WASP (|Pollacco et al.ll2006h . 



Transit searches are most sensitive to large 
planets with short orbital periods, and have now 
found nine Jupiter-mass planets in orbits of less 
than 2 c0. This ha s led to suggestion s of a class of 
very-hot Jupiters ( Melo et aL[2006), with highl y 
irradiated atmospheres (e.g. iFortnev et al.ll2007h . 
to which WASP-4b likely belongs. 

2. Photometric Observations 

The WASP (Wide Angle Search for Planets) 
consortium operates two identical robotic obser- 
vatories; SuperWASP-North on La Palma, in the 
Northern hemisphere, and SuperWASP-South, in 
the Southern hemisphere, hosted by the South 
African Astronomical Observatory (SAAO). Each 
consists of eight wide-field cameras consisting of 
an 11.1-cm aperture Canon 200 mm f/1.8 lens 
backed by a 2kx2k e2v CCD. The eight cameras 
cover 490 square degrees per pointing. WASP-S 
started operating in May 2006, with a strategy of 
tiling six to eight fields with a cadence of 5-10 
mins and exposure times of 30 sees. These fields 
rotate with the seasons, accumulating to a strip 
centered on a declination of -32°, resulting in 
lightcurves for over one million stars brighter than 
13 th magnitude. For details of the WASP project , 



hardw are and data processing see iPollacco et al 
( 20061 ). Further details of our data analysis 



proced ures are given in ICollier Cameron et al 



(|2007aj ). reporting the discovery of WASP-lb 
and WASP-2b from SuperWASP-North, and in 
Pollacco et al] (|2007t ) reporting the discovery of 
WASP-3b. 

From the WASP-S data collected between 2006 
May-November we identified 1SWASP J233415.06- 
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Fig. 1.— The Super WASP-South discovery 
lightcurve of WASP-4 folded on the 1.3-d orbital 
period. 



420341.1 (USNO-B1.0 0479-0948995) as a high- 
priority planetary transit candidate; over 4000 
measurements revealed a possible transit recur- 
ring every 1.3 days (see ICollier Cameron et al. 
20061 for a description of our tr a nsit-se arch 
methods and ICollier Cameron et al.l l2007bl for 
an account of our selection of high-priority 
candidates). A catalogue search within the 48" 
photometric aperture of 1SWASP J233415.06- 
420341.1 revealed no bright blending companions 
or known variable/active stars, which could mimic 
the photometric signature of a transiting planet. 
The WASP-S discovery lightcurve is shown in 
Figure Q] 

A full transit of WASP-4b was observed by 
EulerCAM on the 1.2-m Euler telescope on 2007 
September 25 th . Observations were performed 
in R-band and were heavily de-focussed to allow 
exposure times of ^2 mins, achieving an RMS 
scatter of 1.8 mmags despite poor transparency 
conditions. Part of a transit was also observed 
in SDSS i' band by the 2.0-m Faulkes Telescope 
South (FTS) on 2007 September 27 th . The 
WASP-S, EulerCAM and FTS lightcurves are 
shown in Figure [2l 

3. Spectroscopic Observations with 
CORALIE 

Spectroscopic measurements were obtained 
using the CORALIE spectrograph installed on the 
Euler telescope. CORALIE was originally a twin 
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Fig. 2.— The WASP-S, FTS and EulerCam 
lightcurves showing the transit of WASP-4b. The 
data are shown folded on the orbital period to- 
gether with the best-fitting model determined 
from a simultaneous MCMC fit of the photomet- 
ric and radial-velocity data. The RMS scatter to 
the model fit of the WASP-S, FTS and Euler- 
Cam lightcurves is 15.3, 2.7 and 1.8 mmags re- 
spectively. Phase 1 for the EulerCam lightcurve 
occurred at HJD=2454368. 59121 UT and for the 
FTS lightcurve at HJD=2454371.26766 UT. 



copy of the ELODIE spectrograph ( Baranne et al 
1996). However, triggered by the interest to carry 
out spectroscopic follow-up on transit candidates 
with this instrument, in June 2007 we carried 
out major changes to CORALIE to increase its 
performance on fainter stars. The fibre link and 
the cross-disperser optics have been removed and 
replaced by a new design (Queloz et al. in prep.). 
The double scrambler has also been removed and 
the grism and prism cross-disperser component 
replaced by a series of 4 Schott F2 prisms of 32° 
angle each. The net outcome of this new design 
is to maintain the spectral range from 381 to 
681 nm but with a large efficiency gain of about 
a factor of 6 (8 below 420 nm) and a spectro- 



scopic resolution of 55 000-60000 (increased by 
10-20%). The overall instrumental precision is 
also improved. All exposures are reduced by the 
automated pipeline adapted to the new optical 
configuration. WASP-4 was observed from 2007 
September to November and the radial velocity 
and the vsini computed using a G2-spectral 
template. Exposures were 30 minutes in length, 
achieving a signal-to-noise ratio of ~15. Radial 
velocity variations of semi-amplitude 251 ms _1 
were detected consistent with a planetary-mass 
companion whose orbital period closely matches 
that from the WASP-S transit detections. The 
RV measurements are listed in Table [TJ and are 
shown phase folded and over plotted with the 
best-fitting orbital model in Figure H The RMS 
of the RV measurements to the model fit is 15.3 
ms _1 . 

Table 1: CORALIE RV measurements for WASP- 
4. 

BJD - 2450000 RV (kms^ 1 ) a RV (kms^ 1 ) 



4359.71082 


57.56557 


0.01864 


4362.63121 


57.79551 


0.01883 


4364.65260 


57.67271 


0.02223 


4365.73690 


57.95703 


0.01613 


4372.75799 


57.58619 


0.01496 


4376.68882 


57.64642 


0.01451 


4378.66887 


57.79207 


0.01325 


4379.73630 


57.50846 


0.01470 


4380.61034 


57.77411 


0.01303 


4382.79025 


57.86396 


0.01823 


4383.55277 


57.51000 


0.01434 


4387.61904 


57.49524 


0.01457 


4408.66110 


57.77556 


0.01682 


4409.51932 


57.82666 


0.02156 



3.1. Line-bisector analysis 

Although the amplitude of the RV variation of 
this system is consistent with a planetary com- 
panion, the signal could be mimicked by spectral 
line distortio ns caused by a blended eclipsing 
binary (e.g. ISantos et al.l 120021 ). To confirm the 
presence of a transiting planet it is necessary to 
exclude this scenario. Blended eclipsing binaries 
can be identified from multi-colour photometry 
or through lightcurve modeling ([Torres et al 
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Fig. 3. — CORALIE radial-velocity measure- 
ments folded on the orbital period together with a 
radial-velocity model which includes the expected 
Rossiter-McLaughlin effect for a star with vsini 
= 2 krns" 1 . The RMS of the RV measurements 
to the model fit is 15.3 ms _1 . 



2004), however, the simplest and most reliable 
method is to search for variations in the spectral 
line profiles themselves. This has the added 
advantage of highlighting atmospheric distor- 
tion effects which can also affect the measured 
radial velocities. If the detected radial velocity 
variations are the result of contamination of the 
spectrum by an eclipsing binary, we would expect 
to see distortions in the line profiles in phase 
with the photometric period. The CORALIE 
cross-correlation functions were analysed using 



the line-bisector technique (jQueloz et al.l [2001) 



and no evidence for variation in the bisector spans 
was found, confirming the planetary nature of 
this system. A plot of the bisector spans is shown 
in Figured! 

4. Stellar parameters 

The individual CORALIE spectra have a 
low signal-to-noise ratio, but when co-added 
(and re-binned) give a S/N of around 30-40 
which is suitable for a preliminary photospheric 
analysis of WASP-4. The analysis was performed 
using the UCLSYN spectral synthesis package and 
ATLAS9 models with out convective overshooting 
(ICastelli et al.l Il997l ). The Ha, Na I D and 



Mg I b lines were used as diagnostics of both 
T c ff and logg. The metallicity was estimated 
using the photospheric lines in the 6000-6200 A 



Fig. 4. — Line bisector spans for the CORALIE 
radial velocity measurements showing no signifi- 
cant variation. 

region. However, the co-added spectrum is not of 
sufficient quality to perform a detailed abundance 
analysis. The parameters obtained from this 
analysis are listed in Table [U In addition to 
the spectral analysis, we have also used Tycho 
B, V and 2MASS magnitudes to estimate the 
effective temperature using the Infrared Flux 
Method (jBlackwell fc Shallid Il977h . This gives 
T e g = 5410 ± 240 K, which is in close agreement 
with that obtained from the spectroscopic analy- 
sis. These results imply a spectral type of around 
G7. 

In our spectra the Li I 6708 A line is not 
detected, allowing us to derive an upper limit on 
the Lithium abundance of log n(Li/H) + 12 < 1.0, 
which is slightly less than the Solar photospheric 
value. This implies a minimum age of around 
2 Gyr (jSestito fc Randicbl l2005h . Comparison 
of the temperature and logg wit h the stellar 



Girardi et al 



evolution models of 
maximum-likelihood values M„ 
M Q and R* = 1.15 ± 0.28 R Q . 



(|2000h gives 
= 0.90 ± 0.07 
The distance of 
WASP-4 was calculated as 300 ± 50 pc using the 
angular diameter from the Infrared Flux Method 
and the value of the radius of the star from 
Tabled 

5. Planetary parameters 

The photometric and orbital parameters for 
WASP-4b were determined from a simultaneous 
Markov-Chain Monte-Carlo fit of both the pho- 
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Tabic 2: Stellar parameters of WASP-4 derived 
from a spectral analysis of the CORALIE data. 



Parameter 



Value 



R.A. (J2000.0) 23 h 34 m 15.06 s 
Dec. (J2000.0) -42°03'41.l" 



Teff (K) 

log 9 
[M/H] 
logn(Li) 
usini (kms -1 ) 

Mass (M ) 
Radius (R©) 
Spectral Type 
Distance (pc) 



5500 ± 150 
4.3 ± 0.2 
0.0 ± 0.2 
< 1.0 
2.2 



0.6 
1.0 



0.90 ± 0.07 
1.15 ± 0.28 
G7V 
300 ± 50 



tometric and spectroscopic d ata. The MCMC 



method is described in detail in Collier Cameron et al 



(|2007bf ) and is the same used for the fitting of the 
parameters of WASP-3b (|Pollacco et alJl2007t ). 

The transit lightcurve was modelled using the 
small- planet approximation of Mandel k. Agoll 
(|2002l ). An initial fit of the radial velocity 
measurements showed no evidence for signifi- 
cant orbital eccentricity, as expected for such a 
short-period system. Accordingly the eccentricity 
was set to zero. This leaves seven parameters 
defining the system: epoch (To), orbital period 
(P), duration of the transit (tx) from first to last 
contact, the squared ratio of the planet to stellar 
radius (AF = (Rp/R-k) 2 ), the impact parameter 
(b = a cos i/Rk), the radial velocity amplitude 
(Ki) and the stellar mass M*. Gaussian priors 
were imposed, as determined from the spectro- 
scopic analysis, such that M* = 0.90 ± 0.07 and 
log g = 4.3 ± 0.2. To balance the weight of 
the radial velocities with the photometry in the 
MCMC analysis we added 7 ms _1 of systematic 
error to the errors listed in Table [2 thus obtaining 
a reduced Chi-squared of 1. 

The resulting fitted parameters are listed, 
together with their 1-er uncertainties, in Table [3] 
The parameters derived for the host star are 
consistent with the spectral analysis from the 
CORALIE data. The precision of the transit 



photometry is sufficient to constrain the stellar 
radius and impact parameter such that the uncer- 
tainties in the resultant planetary parameters are 
dominated by the uncertainty in the stellar mass. 

Table 3: Fitted system parameters of WASP- 
4b from a simultaneous MCMC analysis of the 
WASP-S, FTS and EulerCam lightcurves together 
with the CORALIE RV data, assuming a circular 
orbit. 



Parameter 



Value 



Error 



Period (days) 
Epoch (HJD) 
Duration (days) 
(Rp/R*) 2 
b 

i (degrees) 

Ki (kins" 1 ) 
7 (kms -1 ) 
a (AU) 
log g (cgs) 

R* (P©) 
M* (Mq) 

P* Oo) 

-Rp/Rjup 
Mp/M Jup 

PP (PJup) 

log g P (cgs) 
T P (A=0)(K) 



1.3382282 

2454365.91464 

0.0928 

0.0241 

0.13 

88.59 

0.24 

57.7326 

0.0230 

4.45 

0.9370 

0.8997 

1.094 

1.416 
1.215 
0.428 
3.142 
1761 



+0.000003 

-0.000003 

+0.00025 

-0.00023 

+0.0009 

-0.0007 

+0.0005 

-0.0002 

+0.13 

-0.12 

+ 1.36 

-1.50 



+0.01 

-0.01 

+0.002 

-0.001 

+0.001 

-0.001 

+0.016 

-0.029 

+0.04 

-0.03 

+0.077 

-0.072 

+0.038 

-0.085 



+0.068 

-0.043 

+0.087 

-0.079 

+0.032 

-0.044 

+0.023 

-0.034 

+24 

-9 



6. Discussion 

A simultaneous fit of precision photometry and 
radial velocity measurements result in a planetary 
radius of 1.42 Rj up and mass of 1.22 Mj up for 
WASP-4b. This is the second largest transiting 
planet discovered to date , second only to TrES-4b 
(jMandushev et al.l l2007t for plots comparing 



W ASP-4b with other tr ansiting extrasolar planets 
see 



Pollacco et al 2007) 



WASP-4b, with an orbital period of 1.3 d, can 
be compared to other ver y-short-period planet s 
such as 
TrES-3b 



OGLE-TR-56b feonacki et alJ l2003h . 



(O 'Donovan et al. 



20071 ) and WASP-3b 



5 



(jPollacco et al.l I2007T ) . The low rotation velocity 
of the host star WASP-4 is comparable with 
that of TrES-3; in contrast, OGLE-TR-56 and 
WASP-3 have much higher rotation velocities, 
which is thought to be due to the young age of 
these stars. The stars WASP-4 and TrES-3 are 
thought to be older, however the synchronisation 
time for WASP-4 wo uld be longer still at 8 Gyr 
(using the method of lMarcv et al.l[l997l ). so their 
low rotational velocity is consistent with their not 
having been spun-up by their planets. 

A preliminary estimate of the stellar param- 
eters indicates that the parent star is spectral 
type G7V with solar metallicity. Together with 
the short orbital period of 1.34 d this results 
in a blackbody planetary surface temperature, 
assuming zero albedo and isotropic re-radiation, 
of 1760 K. The orbital distance of 0.023 AU places 
WASP-4b well within the criterion (a < 0.04 
AU) for the propose d new class of pM planets 
( Fortnev et al.l 120071 ). These planets display 
a temperature inversion due to low-pressure 
stellar absorption by gaseous TiO and VO. As a 
result they exhibit an unusually hot stratosphere 
which emits strongly in the mid-infrared. The 
WASP-4b system has a larger flux ratio than the 
very-hot-Jupiters WASP-3b and TrES-3b and so 
is an ideal target for secondary eclipse studies 
by Spitzer. The study of more of these systems, 
including WASP-4b, should help constrain models 
for atmospheric dynamics and heat distribution. 
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